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Abstract 



■ In this paper, we show that the expansion history of the Universe in power-law cosmology 
I essentially depends on two crucial parameters, namely the Hubble constant Hq and deceler- 

■ ation parameter q. We find the constraints on these parameters from the latest H{z) and 
SNe la data. At la level the constraints from H(z) data are obtained as q = — O.lS^Q^g 
and Hq = 68.43^2 so ^™ Mpc~^ while the constraints from the SNe la data read as 
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H{z) and SNe la data, which yields the constraints q = — 0.34^qq5 and Hq = 68.931q'52 
km s~^ Mpc~^. The estimates of Hq are found to be in close agreement with some recent 
probes carried out in the literature. The analysis reveals that the observational data success- 



(N : q = -0.38;o;o5 and Hq = 69.18;o:54 km s"^ Mpc~^ We also perform the joint test using 



t^l^ fully describe the cosmic acceleration within the framework of power-law cosmology. We find 

that the power-law cosmology accommodates well the H{z) and SNe la data. We also test 

^Sj ■ the power-law cosmology using the primordial nucleosynthesis, which yields the constraints 

^ ! q> 0.72 and Hq < 41.49 km Mpc~^. These constraints are found to be inconsistent with 

the ones derived from the H{z) and SNe la data. We carry out the statefinder analysis, and 

Q"^ ■ find that the power-law cosmological models approach the standard ACDM model as g ^ —1. 

. Finally, we conclude that despite having several good features power-law cosmology is not a 

G^ ■ complete package for the cosmological purposes. 

O" 
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^ : 1 Introduction 



Power-law cosmology finds a reasonable place in the literature to address some common problems 
(e.g., age problem, flatness problem, horizon problem etc.) associated with the Standard Cold 
Dark Matter (SCDM) model based on Big Bang theory. In such a cosmology, the cosmological 
evolution is described by the scale factor a{t) oc t", where a is a constant. The viability of the 
model with a > 1 has been explored in a series of articles in different contexts (Lohiya & Sethi 
1999; Batra et al. 1999, 2000; Gehlaut et al. 2002, 2003; Dev et al. 2002, 2008; Sethi et al. 2005a, 
2005b; Zhu et al. 2008). Observational constraints on phantom power-law cosmology are discussed 
by Kaeonikhom, Gumjudpai & Saridakis (2011). The motivation for such an endeavor is followed 
by a number of considerations. For instance, power-law cosmological models with a > 1 do not 
encounter the horizon problem at all (Sethi, Dev & Jain 2005). These models do not witness the 
flatness problem since the matter density is not constrained by the scale factor. In these models, 
age of the Universe turns out at least fifty percent greater than the age predicted by the SCDM 
model. This bridges the gap between the age of the Universe and the age estimates of globular 
clusters and high-z redshift galaxies, and thus the age problem is alleviated. Sethi, Dev & Jain 
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(2005) showed that an open hnear coasting cosmological model constrained with type la supernovae 
(SNe la) gold sample and ages of old quasars accommodates a very old high-rcdshift quasar, which 
the SCDM model fails to do. The linear coasting cosmology is found to be consistent with the 
gravitational lensing statistics (Dev et al. 2002) and the primordial nucleosynthesis (Batra et al. 
1999, 2000). Kaplinghat et al. (1999) found that the power-law cosmological models which succeed 
in primordial nuclcosythesis are in conflict with the constraints from Hubble expansion rates and 
SNe la magnitude rcdshift relations. 

In any cosmological model, the Hubble constant Hq and deceleration parameter q play an 
important role is describing the nature of evolution of the Universe. The former one tells us the 
expansion rate of the Universe today while the latter one characterizes the accelerating (q < 0) 
or decelerating (q > 0) nature of the Universe. In the recent past, there have been numerous 
attempts to estimate the value of Hq. Freedman et al. (2001) used the Hubble Space Telescope 
(HST) observations of Ccphcid variable to estimate a value of Hq = 72 ± 8 km s^^ Mpc^^. An 
observational determination of the Hubble constant obtained by Suyu et al. (2010) based upon 
measurements of gravitational lensing by using the HST yielded a value of Hq = 69.7^^0 km s-^ 
Mpc~^. WMAP seven-year results gave an estimate of Hq = 71.0 ± 2.5 km s~^ Mpc~^ based on 
WMAP data alone, and an estimate of Hq = lOA^l'l km s~^ Mpc~^ based on WMAP data with 
Gaussian priors corresponding to earlier estimates from other studies (Jarosik et al. 2011). A 
recent estimate of the Hubble constant, which used a new infrared camera on the HST to measure 
the distance and redshift for a collection of astronomical objects, gives a value of Hq — 73.8 ± 2.4 
km s~^ Mpc~^ (Riess et al. 2011). An alternative probe using data from galactic clusters gave a 
value of Hq = 67.0 ± 3.2 km s'^ Mpc^^ (Beutler et al. 2011). 

In this paper, we show that the power-law cosmology essentially depends on the parameters q 
and Hq. We intend to find the observational constraints on the power- law cosmology parameters 
using the recent observational data from H{z) and supernova observations. We also intend to 
test the power-law cosmology with primordial nucleosynthesis. The paper is organized as follows: 
The basic equations of power-law cosmology are introduced in Section 2. Section 3 deals with the 
constraints on the parameters q and Hq from the latest H{z) data while Section 4 is devoted to 
find the constraints using Union2 compilation of 557 SNe la data points. In Section 5, we perform 
the joint test using H{z) and SNe la data. Section 6 deals with the constraints on power-law 
cosmology from primordial nucleosynthesis while Section 7 is devoted to study the statefinders in 
power-law cosmology. In the last Section, we summarize the main results of the paper. 



We study a general class of power-law cosmology described by the dimensionless scale factor 



where Iq is present age of the universe, Oq is the value of a today and a is a dimensionless posi- 
tive parameter. Hereafter the subscript denotes the present-day value of the parameter under 
consideration. 

The deceleration parameter q, which characterizes accelerating {q < 0) or decelerating (g > 0) 
nature of the Universe, reads as 



2 Basic Equations in Power-law Cosmology 




(1) 
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where an over dot denotes derivative with respect to cosmic time t. The positivity of a leads to 
q>-l. 

The cosmic scale factor in terms of the deceleration parameter may be written as 

^ \ 1/(1+9) 



a{t) = ao . (3) 

We observe that q > —1 is the condition for expanding Universe in the power-law cosmological 
model. 

The expansion history of the Universe is described by the Hubble parameter, 

1 

while the present expansion rate of the Universe is given by Hq 



The scale factor a and the redshift z are connected through the relation a — ao{l + z) ^. 
Therefore, the Hubble parameter in terms of the redshift may be expressed as 

H{z) = Ho{l + zY^'^. (5) 

This shows that the expansion history of the Universe in power-law cosmology depends on the 

parameters Hq and q. 

It may be noted that the above model is well motivated in the literature as mentioned in 
the previous section. However the focus has been on one parameter namely a (or q). But here 
we find the observational constraints on both parameters Hq and q by subjecting the power-law 
cosmological model to the latest data from H{z) and SNe la observations. We also constrain Hq 
and q using primordial nucleosynthesis scenario. 



3 Constraints from observational H{z) data 

Simon, Verde & Jimenez (2005) determined nine H{z) data points in the range < 2; < 1.8 by 
using the differential ages of passively evolving galaxies determined from the Gemini Deep Deep 
Survey and archival data. Recently, H{z) data at 11 different redshifts based on the differential 
ages of red- envelope galaxies were reported by Stern ct al. (2010) while 3 more H{z) data points 
were obtained by Gaztanaga, Cabre & Hui (2009). The newly H{z) data points have been used 
to constrain parameters of various cosmological models (Yang & Zhang 2010; Cao, Zhu & Liang 
2011; Chen & Ratra 2011; Paul, Thakur & Ghose 2010, 2011). Here, we use 13 observational H{z) 
data points given in Table 1 of the paper by Chen & Ratra (2011) and the one at 2; = estimated 
in the work by Riess et al. (2011). For this sake, we define the ^ 
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H{zi, q, Hq) - Hobs{zi 



(6) 



The model has two free parameters namely q and Hq. We perform a grid search in the entire 
parametric space {q > —1 and Hq > 0) to find the best fit model. We find that the best fit values of 
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the parameters are q = —0.18 and Hq = 68.43 together with xl = 1-49, where xl = Xmm/(degree 
of freedom). Here and in what follows Hq is in the units of km s~^ Mpc~^. The negative value 
of q suggests that the power-law cosmological model fitted with the newly obtained H{z) data 
confirms the accelerating nature of the present-day Universe. The likelihood contours at 68.3% 
(inner contour), 95.4% (middle contour) and 99.73% (outer contour) confidence levels around the 
best fit values point (—0.18, 68.43) (represented by star symbol) in the q — Hq plane are shown in 
FigHJ The errors at la level are obtained as g = — 0.18to;i2 — 68-43]^|gQ. 




-0.4 -0.2 0.0 0.2 



Figure 1: The likelihood contours at 68.3% (inner contour), 95.4% (middle contour) and 99.73% (outer contour) 
confidence levels around the best fit values point (—0.18, 68.43) (shown by star symbol) in the q — Hq plane obtained 
by fitting power-law cosmological model with H{z) data. 



4 Constraints from observational SNe la data 

The observations directly measure the apparent magnitude m of a supernova and its redshift z. The 
apparent magnitude m of the supernova is related to the luminosity distance di of the supernova 
through 

m = M + 51og,o(^)+25, (7) 

where M is the absolute magnitude, which is believed to be constant for all SNe la. 
It is convenient to work with Hubble free luminosity distance given by 

Dl{z) = ^dL{z). (8) 

Now, Eq.([7]) can be written as 

m = M + 51ogioI^L(;2)-51ogioi^o + 52.38 . (9) 
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The distance modulus is given by 



fx[z) = m-M = 5 logio Dl{z) - 5 log^Q Hq + 52.38 



(10) 



The Hubble free luminosity distance D^, in the present case, for a geometrically fiat Universe 
reads as 



+ 



Hiz') 



dz' 



1 



[(1 + ^) 
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SNe la are always used as standard candles, and are believed to provide strongest constraints 
on the cosmological parameters. In the present analysis, we use recently released Union2 set of 
557 SNe la from Supernova Cosmology Project (Amanullah et al. (2010). In this case, we define 
the as 



557 



fi{zi,q,Ho) - fiobs{zi 



(J.: 



(12) 



After performing a grid search in the entire parametric space (g > — 1 and > 0), we find that 
the best fit values of the parameters are q = —0.38 and Ho = 69.18 together with xt = 0.99. Again, 
the power-law cosmological model fitted with the 557 SNe la data confirms the cosmic acceleration 
with q = —0.38. The likelihood contours at 68.3% (inner contour), 95.4% (middle contour) and 
99.73% (outer contour) confidence levels around the best fit values point (—0.38,69.18) (indicated 
by star symbol) in the q — Hq plane are shown in FigJ2J The errors at la level are derived as 
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Figure 2: The likelihood contours at 68.3% (inner contour), 95.4% (middle contour) and 99.73% (outer contour) 
confidence levels around the best fit values point (—0.38,69.18) (indicated by star symbol) in the q — Hq plane 
obtained by fitting power-law cosmological model with SNe la data. 
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5 Best fit model from the joint test: i/(z)+SNe la data 



In order to obtain tighter constraints on the model parameters and to avoid degeneracy in the 
observational data, we combine H{z) and SNe la data. Since H{z) and SNe la data are obtained 
from independent cosmological probes, the total likehhood is considered to be the product of 
separate likelihoods of the two probes. Therefore, we define 

xLai =Xh + xIn- (13) 

In the joint analysis, we find that the best fit values of the parameters are q = —0.34 and 
Ho = 68.93 together with xl = 1-01. The joint test also confirms the cosmic acceleration with 
q = —0.34. The likelihood contours at 68.3% (inner contour) and 95.4% (outer contour) confidence 
levels around the best fit values point (—0.34, 68.93) (indicated by star symbol) in the q — Ho plane 
are shown in Fig|31 The errors at Icr level read as g = — 0.341q q5 and Hq = 68.93lo'52- 
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Figure 3: The likelihood contours at 68.3% (inner contour), 95.4% (middle contour) and 99.73% (outer contour) 
confidence levels around the best fit values point (—0.34, 68.93) (shown by star symbol) in the q — Hq plane obtained 
by fitting power-law cosmological model with H(z)+SNe la data. 

FigJH demonstrates the comparison of the best fit cosmological model obtained from the joint 
test with the observational H{z) data in the la region 68.41 < Hq < 69.46, —0.39 < q < —0.29. 
We observe that the model fits well to the observational 15 H{z) points shown with error bars, 
especially at redshifts z < 1. 

In Figini the comparison of the derived best fit model based on iJ(2;)+SNe la with the obser- 
vational 557 SNe la data points (shown with error bars) of Union2 compilation in the la region 
68.41 < Ho < 69.46, -0.39 < q < -0.29 is illustrated. We see that the model is in excellent 
agreement with the observational SNe la data. 
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Figure 4: The observational 15 ll(z) data points are shown with error bars (red color online). Variation of best fit 
model li(z) curve (solid) based on iJ(z)+SNe la data is shown vs z. The dashed curve corresponds to the maximum 
values of ii{z) in the Ict region 68.41 < < 69.46, —0.39 < q < —0.29 while the dotted curve corresponds to the 
minimum values of H{z) in the same region. We observe that the best fit model fits well to the observational data 
points of H{z) especially at redshits z < 1. 
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Figure 5: The observational 557 SNe la data points are shown with error bars (red color online). Best fit 
model distance modulus /i(z) curve (solid) based on H{z)+SNe la data is shown vs z. The dashed curve and the 
dotted curve respectively correspond to the maximum and minimum values of iJ,{z) in the la region 68.41 < Hq < 
69.46, —0.39 < q < —0.29. We observe that the best fit model is in excellent agreement with the observational data 
points of SNe la at all redshifts. 
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6 Constraints from primordial nucleosynthesis 



Before we find constraints on power-law cosmology parameters from primordial big bang nucle- 
osynthesis (BBN) , it is helpful to reproduce the brief review of the BBN in standard model given 
by Kaphnghat et al. (1999). The neutron-proton ratio n/p — exp(— Q/T), where Q = 1.29 MeV 
is the neutron-proton mass difference, at high temperatures T > 1 MeV is maintained by charged- 
current weak interactions among neutrons, protons, electrons, positrons and neutrinos. When the 
Universe is of order 1 s old, T < 1 MeV, the n/p ratio "freezes out" due to the inequilibrium of 
weak interactions and free neutron decay with a lifetime of 887 s. At this stage the deuterium (D) 
produced due to collision of neutrons and protons is rapidly photodissociated by the cosmic back- 
ground photons. This causes very low abundance of D and thus heavier nuclei are not formed at 
this epoch. Thus nucleosynthesis is delayed by this "photodissociation bottleneck" . However, when 
the universe is ~ 3 minutes old, the temperature falls below ~ 80 keV, the deuterium bottleneck 
is broken. At this stage the nuclear reactions quickly burn out the remaining free neutrons into 
^He and leave trace amounts of D, ^He and ^Li (Walker et al. 1991). A viable cosmological model 
must mimic the above scenario for proper synthesis of the light elements in the early Universe. In 
the following, we test the power-law cosmological model for primordial nucleosynthesis. 

In power- law cosmology the scale factor a{t) and the cosmic microwave background temperature 
T{t) are related through the relation: 



where ^ stands for any non-adiabatic expansion due to entropy production. In standard cosmology, 
the instantaneous annihilation is assumed at T = 111^. The heating due to this annihilation is 
accounted by /3 where /3 = 1 for T < rrie while (3 = (11/4)-*^/^ for T > mg. 

In order to find constraints from primordial nucleosythesis on power-law cosmology parameters, 
we utilize ~ 13.7 Gyr (age of the Universe) estimated by Komatsu et al. (2011) on the basis of 7 
year data from WMAP and astrophysical data from other sources. Further, the model is assumed 
to have current temperature Tq = 2.728 K. The primordial nucleosynthesis requires that t < 887 
s when T 80 keV. This puts the following constraints on q and H^: 



We observe that primordial nucleosynthesis requires decelerating expansion (g > 0) of the 
Universe within the framework of power-law cosmology. This is contrary to outcome of accelerated 
expansion (g < 0) of the Universe obtained in earlier sections by using latest data sets of H{z) and 
SNe la observations. Also the primordial nucleosynthesis forces the Hubble constant Hq < 41.49, 
which is much smaller than its values estimated in earlier sections by using the updated data 
sets of H{z) and SNe la observations. This shows that the power- law cosmological models which 
succeed in mimicking the nucleosynthesis scenario are in conflict with the constraints obtained on 
power-law cosmology parameters by using the latest observational data from H{z) and SNe la. A 
similar conclusion on power-law cosmologies was drawn by Kaplinghat et al. (1999). However, in 
the case at hand, it has been done with updated observational data sets and the constraints have 
been obtained on both the parameters q and Hq. 




(14) 



q > 0.72 



and 



Ho < 41.49. 
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7 Constraints on statefinders in power-law cosmology 



The Hubble parameter H = a/a and deceleration parameter q = —d/aH'^ are useful geometric 
parameters in cosmology, which describe the expansion history of the Universe. For instance, 
H > {a > 0) indicates an expanding Universe while g < (a > 0) characterizes an accelerated 
expansion of the Universe. In order to explain accelerated expansion of the Universe, various dark 
energy models have been proposed in the literature. However, these models encounter degeneracy 
on the geometric parameters H (involving first time derivative of scale factor) and q (involving 
second time derivative of scale factor) at the present epoch (Malekjani & Khadom-Mohamaddi 
2012). Thus the geometric parameters H and q are not capable of discriminating between different 
dark energy models and a viable diagnostic tool is required for this purpose. Sahni et al. (2003) 
proposed a pair of parameters {r, s} called statefinders as a means of distinguishing between 
different dark energy models. The statefinders involve derivatives of scale factor upto third order 
and are defined as 

'a r — 1 

and 



aH^ 3{q - 1/2) 

The remarkable feature of statefinders is that these parameters depend on scale factor and its 
time derivatives, and hence are geometric in nature (Sahni 2002). Further, different dark energy 
models exhibit different evolutionary trajectories in the s — r plane. Moreover, the well known flat 
ACDM model corresponds to the point s = and r = 1 in the s — r plane (Alam et al. 2003). 
These features of statefinders provide an opportunity to distinguish between different dark energy 
models. The statefinder diagnostic tool has been extensively used in the literature as a means to 
distinguish between different dark energy models (see Alam et al. 2003; Ali et al. 2010; Malekjani 
& Khadom-Mohamaddi 2012 and references therein). In what follows, we test the power- law 
cosmology by statefinder diagnostic tool. 

In power-law cosmology, the statefinders read as 

r = 2g^ + q and s = |(g + 1), where q ^ 

We immediately notice that r = 1 and s = at g = — 1. Thus the the power-law cosmology 
mimics the ACDM model at g = — 1. We obtain the evolutionary s — r and q — r trajectories for the 
power-law cosmology as shown in Fig. [6]for the values of q in the range —1 < q < 0.5 . The black 
dot in the left panel (a) at (s,r) = (0, 1) represents the location of flat ACDM model while the 
the black dot in the right panel (b) shows the location of the de Sitter (dS) point (g,r) = (—1, 1). 
The star symbol in both panels corresponds to the best fit model based on if(z)+SNe la. It is 
interesting to note that the ACDM statefinder pair (0, 1) or equivalently the dS point (—1, 1) is 
an attractor in power-law cosmology. 

Now we find the constraints on the statefinders of power-law cosmology from H{z), SNe la and 
BBN observations. The H(z) data constrain the statefinders as r = — 0.09]'^o;o3 and s = 0.58]|!o;'j'2 
while the SNe la data constraints on statefinders are r = — 0.09lo;o2 and s = 0.4llg;o3. The joint 
test of H{z) and SNe la data puts the following constraints on statefinders: r = — 0.1llg;o^ and 
s = 0.441q q3. The errors in the above values are at la level. The primordial nucleosynthesis 
restricts r > 1.76 and s > 1.15. We observe that the best fit values of statefinders r and s 
predicted by observational H{z) and SNe la data are in conflict with the ones estimated by BBN, 
as expected. 
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Figure 6: (a) Variation of r versus s. The black dot represents the statefinder pair (s,r) — (0, 1) or equivalently 
location of flat ACDM model in the s — r plane, and star symbol on the s — r curve shows the position of statefinder 
parameters in the best fit model based on iJ(2;)+SNe la data. 

(b) Variation of r versus q. The black dot represents the location of de Sitter (dS) point {q^r) = (—1, 1) in the 
q — r plane, and star symbol on the g — r curve shows the position of q and r in the best fit model obtained from 
H{z)+me la data. 

In both panels, the vertical dashed line separates the acceleration and deceleration zones. The arrows show the 
direction of the evolution of trajectories as q varies from 0.5 to —1. We observe that ACDM point (0, 1) or 
equivalently the dS point (—1,1) is an attractor in power-law cosmology. 
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8 Summary 



In this paper, we have found the bounds on the parameters Hq and q of the power-law cosmology. 
The numerical results are summarized in Table 1. 



Table 1: Summary of the numerical results 



Data/Source 


Q 


Ho 




r 


s 


H{z) 






1.49 


-0.091°:°^ 


0.581H^ 


SNe la 


-0.38tHi 


69.18t°:tt 


0.99 




o.4it°:°i 


ii'(^)+SNe la 


-0.341H^ 




1.01 


-o.iiiS:°? 


44+0-03 


BBN 


>0.72 


< 41.49 




>i.r6 


>1-15 



Some key observations are as follows: 

• The constraints on the deceleration parameter q clearly indicate that the astronomical obser- 
vations of H{z) and SNe la predict the cosmic acceleration within the framework of power- 
law cosmology. Thus, power-law cosmological models are viable for describing the observed 
accelerating nature of Universe. 

• We see that the estimates of Hubble constant in power-law cosmology are in close agreement 
with independent investigations of Hq carried out in literature (Frccdman et al. 2001; Suyu 
et al. 2010; Jarosik et al. 2010; Riess et al. 2011; Beutler et al. 2011) as discussed in Section 
1. 

• The derived best fit model fits well to the observational data points from H{z) and SNe la 
observations (see. Fig. 4 and Fig. 5). 

• The primordial nucleosynthesis demands a decelerating expansion of the Universe with 
smaller values of Hubble constant. 

• The statefinder analysis shows that power-law cosmological models approach the standard 
ACDM model in future (see. Fig. 6) with varying values of q. 

We see that the power-law cosmology turns out viable in the description of the acceleration 
of present-day Universe when subjected to recent observations of H{z) and SNe la. Also the 
Hubble constant within the framework of power-law cosmology falls in the range of observations. 
However, the power-law cosmology fails to produce primordial nucleosytheis with the values of q 
and Hq estimated from observational data of H{z) and SNe la as discussed in Section 6. Moreover, 
because of the constant value of deceleration parameter q in power-law cosmology, it fails to provide 
time or redshift based transition of the Universe from deceleration to acceleration. Thus one has 
to use different values of q for the description of Universe at different epochs. Finally, despite 
having several useful features, the power-law cosmology is not a complete package for cosmological 
purposes. 
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